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Abstract 

We review the present theoretical status of the CP violating ra- 
tio e'/e in the Standard Model (SM) and confront its estimates with 
the most recent data. In particular we review the present status of 
the most important parameters m s , B^ 2 \ B^ 2 \ A^ and Q v + V '- 
While the sign and the order of magnitude of SM estimates for e'/e 
agree with the data, for central values of these parameters most 
estimates of e'/e in the SM are substantially below the grand ex- 
perimental average (21.2 ± 4.6) • 10 -4 . Only in a small corner of 
the parameter space can e'/e in SM be made consistent with exper- 
imental results. In view of very large theoretical uncertainties, it is 
impossible to conclude at present that the data on e'/e indicate new 
physics. A brief discussion of e'/e beyond the Standard Model is 
presented. 
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1 Introduction 



The purpose of this talk is to summarize the present theoretical status of 
the CP violating ratio e'/ e and to confront it with the recent experimental 
findings. The parameters e and e' describe two types of CP violation in 
the decays Kl — > vr7r, which could not take place if CP was conserved. 
In the Standard Model a non-vanishing value of e originates in the fact 
that the mass eigenstate K L is not a CP eigenstates due to the complex 
CKM couplings [[IJ in the box diagrams responsible for the K° — K° mixing. 
Indeed, Kl is a linear combination of CP eigenstates K 2 (CP=— ) and K\ 
( CP=+): Kl = K 2 + eKi, where e is a small parameter. The decay of K L 
via Ki into CP=+ state nn is termed indirectly CP violating as it proceeds 
not via explicit breaking of the CP symmetry in the decay itself but via the 
admixture of the CP state with opposite CP parity to the dominant one. 
The measure for this indirect CP violation is defined as 



_ A(K L -> (7T7r) J=0 ) 
£ A(K S - (7T7r) J= o) : 

which can be rewritten as 



where M 12 represents K° — K° mixing, AMk is the Kl — Ks mass difference 
and A Q is the isospin amplitude for irir in the 1 = state. The second term 
involving £ cancells the phase convention dependence of the first term but 
in the usual CKM convention can be safely neglected. In this limit e = e. 

On the other hand, so-called direct CP violation is realized via a direct 
transition of a CP odd to a CP even state or vice versa. K 2 — > rnr in the 
case of Kl — > vr7r. A measure of such a direct CP violation in Kl — > vr7r is 
characterized by a complex parameter e' defined as 

a* = -Llm (^) exp(i^), $ e , = | + 5 2 - 5 • (1.3) 

Here the subscript 7 = 0, 2 denotes 7T7t states with isospin 0, 2 equivalent 
to AI = 1/2 and AJ = 3/2 transitions, respectively, and <5o,2 are the 
corresponding strong phases. The weak CKM phases are contained in A 
and A 2 . The isospin amplitudes Aj are complex quantities which depend on 
phase conventions. On the other hand, e' measures the difference between 
the phases of A 2 and A Q and is a physical quantity. The strong phases S ^ 2 
can be extracted from -kit scattering. Then $ e / pa tt/4. 
Experimentally we have 0, [!| 

£ cxp = (2.280 ±0.013) ■ 10~ 3 expi$ £ , $ e w (1.4) 
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Until recently the experimental situation on e'/e was rather unclear: 



Re(a'/e)-l (23 ± 7) ' W ^ (NA31) [4] 

ne^/e; | ( 7 4 ± 5 9 ) . 10 -4 (E731) g. 1 J 

While the result of the NA31 collaboration at CERN [|| clearly indicated 
direct CP violation, the value of E731 at Fermilab ||, was compatible with 
superweak theories |J in which e'/e = 0. After heroic efforts on both 
sides of the Atlantic, this controversy is now settled with the two new 
measurements by KTeV at Fermilab f| and NA48 at CERN [g 

r P (V' /f\ - J ( 28 -° ± 41 ) • 10-4 ^ KTeV ) M n 

l / J \ (18.5 ±7.3) ■ 10" 4 (NA48) §, 1 J 

which together with the NA31 result confidently establish direct CP viola- 
tion in nature ruling out superweak models . The grand average including 
NA31, E731, KTeV and NA48 results reads g 

Re(e'/e) = (21.2 ± 4.6) • 1(T 4 (1.7) 

very close to the NA31 result but with a smaller error. The error should 
be further reduced once complete data from both collaborations have been 
analyzed. It is also of great interest to see what value for e' / e will be mea- 
sured by KLOE at Frascati, which uses a different experimental technique 
than KTeV and NA48. 



2 e in the Standard Model 

In the Standard Model e is found by calculating the box diagrams with 
internal u, c, t, W ± exchanges and including short distance QCD correc- 
tions. The final result can be written as follows 

e = C £ B K lm\ t {ReA c faiS^c) - V3 S o( x c%t)} ~ ReX t T) 2 S (x t )} exp(i7r / A) , 

(2.1) 

where C £ is a numerical factor, Aj = V* s Vid and S (xi) with Xi = mf/M^j 
are Inami-Lim functions 0. Explicit expressions can be found in ||10||. The 



NLO values of the QCD factors rji are given as follows [TT 



771 = 1.38 ±0.20, 772 = 0.57 ±0.01, m = 0.47 ± 0.04 . (2.2) 



The main theoretical uncertainty in (|2.1|) resides in the non-perturbative 
parameter Bk which parametrizes the relevant hadronic matrix element 
(K®\{sd)v-A(sd)v-A\K°) . Recent reviews of Bk are given in |1^, [13], [TO 



In our numerical analysis presented below we will use Bk = 0.80 ± 0.15 
which is in the ball park of various lattice and large-N estimates. 
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It is well known that that the experimental value of e in ( |1.4| ) can be 
accomodated in the Standard Model. We know this from numerous analyses 
of the unitarity triangle which in addition to ( |1.4[ ) take into account data 
on B^ s — B® s mixings and the values of the CKM elements \ V US \, \V u b\ and 
\V c b\. From this analysis one extracts in particular 



f (1.33 ± 0.14) • 10~ 4 (Monte Carlo) . . 

JmA t - j (1 33 ± Q 3Q) _ 10 _ 4 ( Scanriirig ). \ Li ) 

The "Monte Carlo" error stands for an analysis in which the experimentally 
measured parameters, like \V US \, \V u b\ and \V c b\, are used with Gaussian 

errors and the theoretical input parameters, like B K , F b \Jb~b are taken 
with flat distributions. In the "scanning" estimate all input parameters are 
scanned independently in the appropriate ranges. Details can be found in 



This analysis gives also sin 2(3 = 0.73 ± 0.09 (0.71 ± 0.13), where (3 is 
one of the angles in the unitarity triangle. The recent study of CP violation 
in B° — ► J /ipKs by CDF fl4| gives sin 2(3 = 0.79 ± 0.44 in good agreement 



with the value above, although the large experimental error precludes any 
definite conclusion. 



3 e' Is in the Standard Model 



3.1 Basic Formulae 

In the Standard Model the ratio e'/e is governed by QCD penguins and 
electroweak penguins (7 and Z° penguins). With increasing value of m t the 



electroweak penguins become increasingly important UT5| , |16|| and entering 
e'/e with the opposite sign to QCD penguins, suppress this ratio for large 
m t . The size of e'/e is also strongly affected by QCD renormalization group 
effects. Without these effects e'/e would be at most O(10 -5 ) independently 



of m t Wfi 



The parameter e' is given in terms of the isospin amplitudes Aj in ( |1.3| ) . 
Applying the operator product expansion to these amplitudes one finds 



where 
with 



ImAt ■ F £ ,{m u A§> m s , Bf>, Bf',^, 



(1/2) D (3/2) 



MS 



p(l/2) _p(3/2)l exp(i$), $ = $ e /-$ £ 



P {l/2) = r5>(A*)<Qi(A0>o(l-n 
P (3/2) = L^y^m^h . 



(3.1) 
(3.2) 

(3.3) 
(3.4) 
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Here 



G-pOJ 



2\e\ReAo 



(Qi)j = ((Tr^rlQilK) 



LO 



ReA 2 
ReA 



(3.5) 



and \i is the renormalization scale which is (9(1 GeV). Since $ ~ 0, F e i 
and e'/e are real to an excellent approximation. The arguments of F £ i will 
be discussed shortly. Explicit expressions for the operators Qi,2 (current- 
current), Q3-6 (QCD-penguins) and Q7-10 (electroweak penguins) are given 
in 111 Oil . The dominat are these two: 



a dp)v-A 



a)V+Ai 



-U,d,: 



3 (- 



d 



P)V-A 



E 



e q (q/3Qa)v+A 



(3.6) 

where, e q denotes the electric quark charges. 

The Wilson coefficient functions yi(fi) were calculated including the 

The de- 
Their 



complete next-to-leading order (NLO) corrections in fnj|-|r§ 



tails of these calculations can be found there and in the review 20 



numerical values for corresponding to a- M< 
two renormalization schemes (NDR and HV) can be found in Rl3f| . 

It is customary in phenomenological applications to take ReAo and u 
from experiment, i.e. 



.,(5) 



(M z 



0.119 ± 0.003 and 



ReA = 3.33- 10" ' GeV 



LO 



0.045, 



(3.7) 



where the last relation reflects the so-called AI = 1/2 rule. This strategy 
avoids to a large extent the hadronic uncertainties in the real parts of the 
isospin amplitudes Aj. 

The sum in (|3.3j ) and (|3.4| ) runs over all contributing operators. p( 3 / 2 ) is 
fully dominated by electroweak penguin contributions. P^ 1 / 2 ) on the other 
hand is governed by QCD penguin contributions which are suppressed by 
isospin breaking in the quark masses (m u 7^ mj. The latter effect is 
described by 



n 



77+77' 



1 (1mA 



2 LB. 



UJ 



ImAn 



0.25 ±0.08. 



(3. 



3.2 History of e'/e 

There is a long history of calculations of e'/ e in the Standard Model. As it 



has been already described in [10], [I3|] , I will be very brief here. The first 
calculation of e'/e for m t Mw without the inclusion of renormalization 
group effects can be found in [p4 . Renormalization group effects in the 



leading logarithmic approximation have been first presented in [ 25| . For 
m t M w only QCD penguins play a substantial role. With increasing m t 
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the role of electroweak penguins becomes important. The first leading log 
analyses for arbitrary m t can be found in [15, 16], where a strong cancella- 
tion between QCD penguins and electroweak penguin contributions to e'/e 
for m t > 150 GeV has been found. Finally, during the nineties consider- 
able progress has been made by calculating complete NLO corrections to 

Q-0 



the Wilson coefficients relevant for e (liTl and e' 



The progress in 

calculating the corresponding hadronic matrix elements was substantially 
slower. It will be summarized below. 

Now, the function F e > in (|3.1|) can be written in a crude approximation 
(not to be used for any serious analysis) as follows 



13- 



110 MeV 
m s (2 GeV) 



A. 



MS 



(1/2), 



340 MeV 

(3/2) 



m t 



165 GeV 



2.5 



(3.9) 



Here Bi are hadronic parameters defined through 

Wo = Bi m (Qe)t c) , (Qsh = B$ l2) (Q 8 )t c) 



(3.10) 



The label "vac" stands for the vacuum insertion estimate of the hadronic 
matrix elements in question for which B^ 2 ^ = B^ 2 ^ = 1. The dependence 
on m s in (|3~9| ) originates in the m s dependence of (Qe}^^ an d {Qs)2 &C \ 
so that {Qq)q and {Q&)2 are roughly proportional to 



Rp. = B, 



(1/2) 



137 MeV 



m s (m c ) + m d (m c ) 



Rs = B 



(3/2) 



137 MeV 



m s (m c ) + m A {m c ) 

(3.11) 

respectively. The scale m c = 1.3 GeV turns out to be convenient. m s (m c ) 



1.17 • m s (2GeV). The formula (|3.9|) exhibits very clearly the dominant 



uncertainties in F F > which reside in the values of m s , Bp, 



(1/2) B (3/2) A (4) 



A^and 



Q v +r]'- Moreover, the partial cancellation between QCD penguin (B^ 2 ^) 
and electroweak penguin [B% ) contributions requires accurate values of 



-Bg 1 ^ and B^ /z> for an acceptable estimate of e'/e. Because of the accurate 
value m t (m t ) = 165 ± 5 GeV, the uncertainty in e'/e due to the top quark 
mass amounts only to a few percent. A very accurate analytic formula for 
F e i has been derived in |26| and its update can be found in [[L3| . 

Now, it has been known for some time that for central values of the 
input parameters the size of e'/e in the Standard Model is well below the 
NA31 value of (23.0 ± 6.5) ■ 10~ 4 . Indeed, extensive NLO analyses with 



>(3/2) 



lattice and large-N estimates of 5 6 ~ 1 and B% ' ' 1 performed first 
in [ IS], |19| and after the top discovery in [27, ^8], [29| have found e' / e in the 
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ball park of (3 - 7) • KT 4 for m s (2 GeV) « 130 MeV. On the other hand 
it has been stressed repeatedly in [[[(], |28| that for extreme values of B^ 2 " 1 , 
Bg and m s still consistent with lattice, QCD sum rules and large-N 
estimates as well as sufficiently high values of lm\ t and A^-, a ratio e'/e 
as high as (2 — 3) • 10~ 3 could be obtained within the Standard Model. Yet, 
it has also been admitted that such simultaneously extreme values of all 
input parameters and consequently values of e'/e close to the NA31 result 
are rather improbable in the Standard Model. Different conclusions have 
been reached in ||30|| , where values (1 — 2) • 10 -3 for e'/e can be found. Also 
the Trieste group PH, which calculated the parameters B^ 2 ^ and B^ 2 ^ 



in the chiral quark model, found e'/e = (1.7tJ; 4 ) • 10~ 3 . On the other hand 
using an effective chiral lagrangian approach, the authors in [32] found e'/ e 
consistent with zero. 



3.3 Hadronic Matrix Elements 

The main source of uncertainty in the calculation of e'/e are the hadronic 
matrix elements {Qi(n))i- They generally depend on the renormalization 
scale fi and on the scheme used to renormalize the operators Qi. These 
two dependences are canceled by those present in the Wilson coefficients 
yi(n) so that the resulting physical e'/e does not (in principle) depend 
on /i and on the renormalization scheme of the operators. Unfortunately, 
the accuracy of the present non-perturbative methods used to evalutate 
{Qi)i is not sufficient to have the \i and scheme dependences of {Qi)i fully 
under control. We believe that this situation will change once the non- 
perturbative calculations, in particular lattice calculations improve. 

As pointed out in the contributions of (V — A) <g> (V — A) operators 
(Qi with i=l,2,3,4,9,10) to e'/e can be determined from the leading CP 
conserving K — > ttti decays, for which the experimental data is summarized 
in ( |3.7| ). The details of this approach will not be discussed here. For the 
central value of lm\ t these operators give a negative contribution to e'/e 
of about —2.5 ■ 10 -4 . This shows that these operators are only relevant if 
e'/e is below 1 ■ 10 -3 . Unfortunately the matrix elements of the dominant 
(V — A) <g> (V + A) operators {Qq and Q$) cannot be determined by the CP 
conserving data and one has to use non-perturbative methods to estimate 
them. Let us then briefly review the present status of the corresponding 
non-perturbative parameters B^ 2 ^ and B^ 2 ^ as well as of m s , £l v+r] i and 
A^ which all enter the evaluation of e'/e as seen in (|3.9|) . 

3.4 The Status of m s , B^ 1/2) , B§ /2 \ Vt^ and 
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3.4.1 



nu 



The values for m s (2 GeV) extracted from quenched lattice calculations and 
QCD sum rules before summer 1999 were 



m s (2GeV) 



110 ±20) MeV 
124 ±22) MeV 



(Lattice) [12 
(QCDS) 



(3.12) 



The value for QCD sum rules is an average |H| over the results given in |34 



QCD sum rules also allow to derive lower bounds on the strange quark mass. 
It was found that generally m s (2GeV) ^ 100 MeV [[35]. The most recent 
quenched lattice results: m s (2GeV) 

(2 GeV) 



97 ±4 MeV 37 



106 ± 7 MeV 
105 ± 5 MeV 



[2 GeV) 



38 are consistent with 



these bounds as well as ( |3.12| ) but have a smaller error. The unquenching 
seems to lower these values down to m s (2GeV) = 84 ± 7 MeV [j39fl . We 
refer to the talks of Ryan and Martinelli [fHJ for more details. The most 
recent determination of m s from the hadronic r-spectral function j|l[] reads 
m s (2GeV) = (114 ± 23) MeV || in a very good agreement with (|1|) 
and lower that the value m s (2GeV) = (170tff) MeV obtained by ALEPH 



43| using this method early this year. We conclude that the error on m s 



decreased considerably during the last two years and the central value is 
in the ball park of m s (2 GeV) = 110 MeV with smaller values coming from 
unquenched lattice QCD. 



3.4.2 B, 



(1/2) 



and B 



(3/2) 



We recall that in the large-N limit B, 



(1/2) 



B 



(3/2) 



1 |4j, |45[. The values 



for B\ L '*> and fl^ obtained in various approaches are collected in table [l|. 
The lattice results have been obtained at \i — 2 GeV. The results in the 
large-N approach |^4|, ^5| and the chiral quark model correspond to scales 
below 1 GeV. However, as a detailed numerical analysis in [|11| showed, 

(1/2) (3/2) 

Bq ' ' and Bp"' are only weakly dependent on fi. Consequently the com- 
parison of these parameters obtained in different approaches at different \x 
is meaningful. Next, the values coming from lattice and chiral quark model 
are given in the NDR renormalization scheme. The corresponding values 
in the HV scheme can be found using approximate relations JT3 



1-^6 Jhv 



1.2(5, 



(1/2) > 



WDR, 



(3/2) > 



mv 



l.2(B£ /2) ] 



NDR- 



(3.13) 



The present results in the large-N approach are unfortunately not sensitive 
to the renormalization scheme but this can be improved 0, 

Concerning the lattice results for B^ 2 \ the old results where in the 
range 0.6 — 1.2 H57|. However, a recent work in ||48H shows that lattice cal- 



culations of B^ 2 ^ are very uncertain and one has to conclude that there are 
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no solid predictions for B^ 2 ^ from the lattice at present. The average value 
of Sg 1//2 ' ) in the large-N approach including full p 2 and p°/N contributions 
and given in table p] is close to 1.0 where the uncertainty comes from the 
variation of the cut-off A c in the effective theory. On the other hand, it has 
been found [54], [5l| that a higher order term 0(p 2 /N) enhances B^ 2 ^ to 
1.5. This result is clearly interesting. Yet, in view of the fact that other 
p 2 /N terms as well as p A and p°/N 2 terms have not been calculated, it 
is premature to take this enhancement seriously. Finally, the chiral quark 
model gives in the NDR scheme the value for 5g as high as 1.33 ± 0.25. 



Table 1: Results for Bq and B$ obtained in various approaches. 



Method 


R (l/2) 


R (3/2) 


Lattice H, 51, 50] 
Large-N|53|, 52] 
ChQM0 


0.72 - 1.10 
1.07-1.58 


0.69-1.06 
0.42 - 0.64 
0.75 - 0.79 



The status of B^ 2 ^ looks better. Most non-perturbative approaches 
find B^ 2 ^ below unity. The suppression of B^ 2 ^ below unity is rather 
modest (at most 20%) in the lattice approaches and in the chiral quark 
model. It is stronger in the large-N approach. Interestingly in the latter 

('1/2) (3/2) 

approach the ratio B 6 / B% ~ 1.72 independently of the cut-off A c . 

Guided by the results presented above and biased to some extent by 
the results from the large-N approach and lattice calculations, we will use 
m our numerical analysis below B<£ /2) and Bf /2) in the ranges: 

B^ /2) = 1.0 ± 0.3, Bf /2) = 0.8 ± 0.2 (3.14) 

keeping always B^ 2 ^ > B^ 2 \ 

3.4.3 a n+v , and A§jL 



The last estimates of £l, n+ri i have been done more than ten years ago [£1 
p3fl and it is desirable to update these analyses which can be summarized 
by Q v +r]' = 0.25 ± 0.08 . In the numerical analysis presented below we 
have incorporated the uncertainty in tt v + v > by increasing the error in B^ 2 ^ 
from ±0.2 to ±0.3. Finally from a s (M z ) = 0.119 ± 0.003 one deduces 
Ag^ = (340 ± 50) MeV. 

3.5 Comments 
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3.5.1 General Remark 



We would like to emphasize that it would not be appropriate to fit Bq , 
B^ 2 \ m s , Aj^L, Q v+r] > and Bk in order to make the Standard Model com- 
patible simultaneously with experimental values on e'/e, e and the analysis 
of the unitarity triangle. Such an approach would be against the whole 
philosophy of searching for new physics with the help of loop induced tran- 
sitions as represented by e'/e and e. Moreover it would not give us any clue 
whether the Standard Model is consistent with the data on e'/e. Indeed, 
it should be kept in mind that: 

• B^ 2 \ B^ 2 \ Bk and {l v +ri' in spite of carrying the names of non- 
perturbative parameters, are really not parameters of the Standard 
Model as they can be calculated by means of non-perturbative meth- 
ods in QCD. 

• m s , A^-, m t , I Kb I and \V u b\ are parameters of the Standard Model 
but there are better places than e'/ e to determine them. In particular 
the usual determinations of these parameters can only marginally be 
affected by physics beyond the Standard Model, which is not neces- 
sarily the case for e and e'/e. 

Consequently, the only parameter to be fitted by direct CP violation is 
lm\ t . The numerical analysis of e'/e as a function of B^ 2 \ B^ 2 \ m s and 
A^ presented below, should only give a global picture for which ranges of 
parameters the presence of new physics in e'/e should be expected. 



3.5.2 The Issue of Final State Interactions 



In ( |1.3|) and (|3.2|) the strong phases 8q ~ 37° and 82 ~ —7° are taken from 
experiment. They can also be calculated from NLO chiral perturbation 
for 7T7T scattering pTfl . However, generally non-perturbative approaches to 
hadronic matrix elements are unable to reproduce them at present. As 81 
are factored out in ( |1.3| ), in non-perturbative calculations in which some 
final state interactions are present in (Qi)i one should make the following 
replacements in (OF) and (|3.4|): 



(Qih 



(cos 8 



nth 



(3.15) 



in order to avoid double counting of final state interaction phases. Here 
(cos<5f)th is obtained in a given non-perturbative calculation. Yet, in most 
calculations the phases are substantially smaller than found in experiment 
44], m, H, H and (Qi)j « Re(Qi)j. 
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The above point has been first discussed by the Trieste group |JT] who 
suggested that in models in which at least the real part of (Qi)i can be 
calculated reliably, one should use (cos5/) exp in fl3.15| ). As (cos^ )ex P ~ 0.8 
and (cos#2)exp ~ 1 this modification enhances P^ 1 / 2 ) by 25% leaving p( 3 / 2 ) 
unchanged. To our knowledge there is no method for hadronic matrix 
elements which can provide 5o ~ 37° and consequently this suggestion may 
lead to an overestimate of the matrix elements and of e'/e. 



3.5.3 e'/e and the AI = 1/2 Rule 



In one of the first estimates of e'/e, Gilman and Wise |25[] used the sug- 
gestion of Vainshtein, Zakharov and Shifman |58| that the amplitude KeAo 



is dominated by the QCD-penguin operator Q G . Estimating (Qq)o in this 
manner they predicted a large value of e'/e. Since then it has been under- 
stood HH, ^2| that as long as the scale \i is not much lower than 1 GeV 
the amplitude ReA is dominated by the current-current operators Qi and 
Q2, rather than by Qq. Indeed, at least in the HV scheme the operator 
Qq does not contribute to ReA for /i = m c at all, as its coefficient z G (m c ) 
relevant for this amplitude vanishes. Also in the NDR scheme ZQ(m c ) is 
negligible. 

For decreasing /i the coefficient z$ (/x) increases and the Qq contribution 
to KeA is larger. However, if the analyses in [Ell HR E2] are taken into 



account, the operators Qi and Q 2 are responsible for at least 90% of ReA 
if the scale \i = 1 GeV is considered. Moreover, it should be stressed that 
whereas Qs operator is irrelevant for the AI =1/2 rule, it is important for 
e'/e. Similarly, whereas the Wilson coefficients ye(fj) and t/sCa*) entering e'/e 
can be sensitive to new physics as they receive contributions from very short 
distance scales, this is not the case for Ze(fj), which due to GIM mechanism 
receives contribution only from fi < m c . Therefore, even if (Qe)o enters 
both e'/e and ReAo, there is no strict relation between the large value of 
e'/e and the AI = 1/2 rule as sometimes stated in the literature. On the 
other hand if the long distance dynamics responsible for the enhancement 
°f (Q 1,2)0 also enhances (Qe)o, then some connection between e'/e and the 
AI = 1/2 rule is possible. There are some indications that this indeed 
could be the case [01], 03, 151 



3.5.4 e'/e, B ( 6 1/2) , B\ i/Z) and m s 

At this symposium there has been a vigorous discussion whether the value of 
m s is really relevant for the estimate of e'/ e. In non-perturbative approaches 
in which hadronic matrix elements can only be calculated in terms of m s , 
F n , Fx etc., it is obvious that the value of m s enters the estimate of e'/e. 
This is in particular the case of the large-N approach. Moreover in this 
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approach the values of B^ 2 ^ and B^ 2 ^ are independent of m s 



On the other hand, as seen in (|3.9| ), e'/e depends approximately on 
B^ 2 \ 5g 3//2 ' ) and m s only through Rq and Rs defined in (|3.11|) . If a given 
non-perturbative approach is able to calculate directly the relevant hadronic 
matrix elements and consequently Rq and i?§ then in principle the value 
of m s may not matter. If indeed Rq and Rs are independent of m s then 



through Q3.ll]) there must be a quadratic dependence of B^ 2 ^ and Bg~' 



J 8 

on m s . As stated above this dependence is not observed in the large- 
N approach |44|, |53|, |52fl . As discussed by Martinelli at this symposium, 



this question is being investigated in the lattice approach at present. On 
the other hand there are results in the literature showing a strong m s - 
dependence of the Bi parameters. This is the case for 5g in the chiral 
quark model where B^ 2 ^ scales like m s [[H]. Similarly values for B^ 2 ^ 
calculated in |60| show a strong m s - dependence. However, B^ 2 ^ calculated 



in the chiral quark model shows a very weak dependence on m s . This 
discussion shows that comparision of (Bi 1/2 \ B^ 2) ) obtained m various 
approaches has to be done with care. 

Personally, I think that the value of m s is relevant for the matrix ele- 

(1/2) (3/2') 

ments (Qe/o and (Qs)2 and (£>g ,Bg ) are nearly m s independent. As 
both operators (after Fierz transformation) have the density-density struc- 
ture, their matrix elements are proportional to the square of the quark 
condensate and hence proportional to 1/m 2 . In this context it should be 
recalled that the anomalous dimenions of Qq and Q$ are in a good ap- 
proximation equal twice the anomalous dimension of the mass operator. 
As a result of this, the products Vq{ij){Q&{ij))o and ys(l-i)(Q 's(/-0)2 and the 
corresponding contributions to e'/e are only very weakly //-dependent. 

3.6 Numerical Estimates of e'/e 
3.6.1 Munich Analysis 

We will begin by presenting the analysis in [|l3j. In table |2| we summarize 
the input parameters used there. The value of m s (m c ), corresponds roughly 
to m s (2 GeV) = (110 ±20) MeV as obtained in lattice simulations. ImA t is 
given in ( |2.3j ) except that in evaluating e' / e the correlation in m t between 
ImA t and F £ i has to be taken into account. In what follows we will present 
the results of two types of numerical analyses of e' j e which use Monte Carlo 
and scanning methods discussed already in section 2. 

Using the first method we find the probability density distributions for 
e' j e in fig. |I|. From this distribution we deduce the following results: 

1 (NDR) (3 16) 
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Table 2: Collection of input parameters. m s (m c ) = 1.17m s (2 GeV). 



Quantity 


Central 


Error 


Reference 


MS 

m s (m c ) 

R (l/2) 

d(3/2) 
■°8 


340 MeV 
130 MeV 

1.0 

0.8 


±50 MeV 
±25 MeV 

±0.3 

±0.2 


i 61] 
See Text 
See Text 
See Text 



The difference between these two results indicates the left over renormaliza- 
tion scheme dependence. Since, the resulting probability density distribu- 
tions for e'/e are very asymmetric with very long tails towards large values 
we quote the medians and the 68%(95%) confidence level intervals. Using 
the second method we find : 

1.05 • 10~ 4 < e'/e < 28.8 • 10" 4 (NDR). (3.17) 

and 

0.26 • 10~ 4 < e'/e < 22.0 • 10" 4 (HV). (3.18) 

We observe that e'/e is generally lower in the HV scheme if the same 
values for B^ 2 ^ and B^ 2 ^ are used in both schemes. Since the present 
non-perturbative methods do not have renormalization scheme dependence 

(1/2) (3/2) 

fully under control we think that such a treatment of Bq and B$ is the 
proper way of estimating scheme dependences at present. Assuming, on the 
other hand, that the values in (|3.14|) correspond to the NDR scheme and 
using the relation ( |3.13| ), we find for the HV scheme the range 0.58 • 10 -4 < 
e'/e < 26.9 ■ 10 -4 which is much closer to the NDR result in Q3.17 ). This 



exercise shows that it is very desirable to have the scheme dependence under 
control. 

We observe that the most probable values for e'/e in the NDR scheme 
are in the ball park of 1 • 10~ 3 . They are lower by roughly 30% in the HV 
scheme if the same values for (Bg , B% ) are used. On the other hand 
the ranges in ( p. 17] ) and ( |3.18| ) show that for particular choices of the input 
parameters, values for e'/e as high as (2 — 3) • 10~ 3 cannot be excluded at 
present. Let us study this in more detail. 

In table ^| we show the values of e'/e in units of 10~ 4 for specific val- 
ues of Bi 1/2 \ Bf /2) and m s (m c ) as calculated in the NDR scheme. The 
corresponding values in the HV scheme are lower as discussed above. The 
fourth column shows the results for central values of all remaining parame- 
ters. The comparison of the the fourth and the fifth column demonstrates 
how e'/e is increased when A^ is raised from 340 MeV to 390 MeV. As 
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-5 5 10 15 20 25 30 35 



e'/e [10- 4 ] 

Figure 1: Probability density distributions for e'/e in NDR and HV schemes 



stated in ( |3.9| ) e'/e is roughly proportional to A^. Finally, in the last 
column maximal values of e'/e are given. To this end we have scanned all 
parameters relevant for the analysis of ImAt within one standard deviation 
and have chosen the highest value of A^L = 390 MeV. Comparison of the 
last two columns demonstrates the impact of the increase of ImA t from its 
central to its maximal value and of the variation of writ. 

Table ^| gives a good insight in the dependence of e' / e on various pa- 
rameters which is roughly described by ( |3.9[ ). We observe the following 
hierarchies: 

• The largest uncertainties reside in m s , B^ /2) and Bf m . 

• The combined uncertainty due to ImAt and m t , present both in ImA t 
and F £ i, is approximately ±25%. The uncertainty due to m t alone is 
only ±5%. 



The uncertainty due to is approximately ±16%. 



• The uncertainty due to Q v+V > is approximately ±12%. 

The large sensitivity of e'/e to m s has been known since the analyses 
in the eighties. In the context of the KTeV result this issue has been first 
analyzed in [p2]| . It has been found that provided 2B^ 2 ^ — B^ 2 ^ < 2 
the consistency of the Standard Model with the KTeV result requires the 
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2a bound m s (2GeV) < 110 MeV. Our analysis is compatible with these 
findings. 



Table 3: Values of e'/e in units of 10 4 for specific values of B$ 2 \ Bg 
and m s (m c ) and other parameters as explained in the text [p~3]j . 



D 


a 


m s (m c )\MeV] 

a\ W L J 


Central 


Ag^ = 390 MeV 


Maximal 






105 


20.2 


23.3 


28.8 


1.3 


0.6 


130 


12.8 


14.8 


18.3 






155 


8.5 


9.9 


12.3 






105 


18.1 


20.8 


26.0 


1.3 


0.8 


130 


11.3 


13.1 


16.4 






155 


7.5 


8.7 


10.9 






105 


15.9 


18.3 


23.2 


1.3 


1.0 


130 


9.9 


11.5 


14.5 






155 


6.5 


7.6 


9.6 






105 


13.7 


15.8 


19.7 


1.0 


0.6 


130 


8.4 


9.8 


12.2 






155 


5.4 


6.4 


7.9 






105 


11.5 


13.3 


16.9 


1.0 


0.8 


130 


7.0 


8.1 


10.4 






155 


4.4 


5.2 


6.6 






105 


9.4 


10.9 


14.1 


1.0 


1.0 


130 


5.5 


6.5 


8.5 






155 


3.3 


4.0 


5.2 



3.6.2 Comparision with Other Analyses 

In table [| we compare the results discussed above with the most recent 
results obtained by other groups. Details of these analyses can be found in 
other contributions to this symposium and in references quoted in the table. 
There exists no recent phenomenological analysis from the Trieste group 
31|] and we quote their 1998 result. The values for B^ 2 ^ and B^ 2 ^ are 



given in the NDR scheme. The corresponding values in the HV scheme can 
be found by using ( |3.13| ). All groups use the Wilson coefficients calculated 



m 



17j-||19| and the differences in e'/e result dominantly from different 
values of B^ 2 \ B^ 2 ^ and m s and to some extend in different values in 
input parameters (like B^) needed for the determination of lm\ t . The 
results for e'/e are given in the NDR scheme except for the result from 
Trieste which corresponds to the HV scheme. The values of m s (m c ) used 
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by all groups are very close to the ones given in table ^|. 

We observe that the results from Munich, Rome and Dortmund are com- 
patible with each other and generally well below the experimental data. The 
values of e'/e for central values of parameters obtained by the Dortmund 
group are in good agreement with those of the Munich group. On the other 
hand as discussed above and also found by the Dortmund group for extreme 
choices of input parameters values for e'/e consistent with experiment can 
be obtained. In the framework of an effective chiral lagrangian approach 
H B^ /2) and Bf /2) cannot be calculated. For B^ /2) = Bf ,2) = 1 values 
for e'/e consistent with zero are obtained. 

The Trieste group finds generally higher values of e' / e, with the central 
value around 17 • 10~ 4 and consequently consistent with the experimental 
findings. Similarly the Dortmund group finds the central values for e'/e in 
the ball park of 15 ■ 10 -4 if they use B^ 2 ^ m 1.5 as obtained by including 
one p 2 /N term. 

Table 4: Results for e'/e in units of 10~ 4 obtained by various groups. 
The labels (MC) and (S) in the last column stand for "Monte Carlo" and 
"Scanning" respectively, as discussed in the text. 



Reference 


£(1/2) 


Bf ,2) 


e'/e[\0-*} 


Munich [13| 
Munich PI 


1.0 ±0.3 
1.0 ±0.3 


0.8 ±0.2 
0.8 ±0.2 


7.71^ (MC) 
1.1 -> 28.8 (S) 


Rome |56|| 


0.83 ±0.83 


0.71 ±0.13 


4.7^3 (MC) 


Dortmund J54J 


0.72 - 1.10 


4 1/2) /i-72 


2.1 -> 26.4 (S) 


Trieste 


1.33 ±0.25 


0.77 ±0.02 


7^ 31 (S) 


Dubna-DESY H 


1.0 


1.0 


-3.2 -> 3.3 (S) 



4 e'/e Beyond the Standard Model 

As we have seen the estimates of e'/e within the Standard Model (SM) are 
generally below the data but in view of large theoretical uncertainties stem- 
ming from hadronic matrix elements one cannot firmly conclude that the 
data on e'/e indicate new physics. However there is still a lot of room for 
non-standard contributions to e'/e and the apparent discrepancy between 
the SM estimates and the data invites for speculations. Indeed results 
from NA31, KTeV and NA48 prompted several analyses of e'/e within var- 
ious extensions of the Standard Model like general supersymmetric models 
[ p2] , |B4| , BH j , models with anomalous gauge couplings |6B| and models with 
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additional fermions and gauge bosons [70]. Unfortunately several of these 
extensions have many free parameters and are not very conclusive. 

On the other hand it is clear that the e'/e data puts models in which 
there are new positive contributions to e and negative contibutions to e' in 
serious difficulties. In particular as analyzed in |L3| the two Higgs Doublet 
Model II can either be ruled out with improved hadronic matrix elements 
or a powerful lower bound on tan/3 can be obtained from e'/e. In the 
Minimal Supersymmetric Standard Model, in addition to charged Higgs 
exchanges in loop diagrams, also charginos contribute. For suitable choice 
of the supersymmetric parameters, the chargino contribution can enhance 
e'/e with respect to the Standard Model expectations fl6£H . Yet, gener- 
ally the most conspicuous effect of minimal supersymmetry is a depletion 
of e'/e. The situation can be different in more general models in which 
there are more parameters than in the two Higgs doublet model II and 
in the MSSM, in particular new CP violating phases. As an example, in 
general supersymmetric models e'/e can be made consistent with experi- 
mental findings through the contributions of the chromomagnetic penguins 
H], ^2], |64], |65| and enhanced Z°-penguins with the opposite sign to the one 
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While substantial new physics contributions to e'/e from chromomag- 
netic penguins and modified Z-penguins appear rather plausible, one can 
give rather solid arguments that new physics should have only a minor im- 
pact on the QCD-penguins represented by the operator Qq. The point is 
that the contribution of Qq to e'/e can generally be written as follows 



e'/e) Qe « Rq [11 - 1.3E(m u ..)] (4.19) 



where 



Rq = ImA* 



110 MeV 


2 


( \ 

MS 


_m s (2 GeV)_ 




I 340 MeV j 



B£ /2) (4.20) 



and E(m t , ...) results from QCD penguin diagrams in the full theory which 
in addition to VT-boson and top-quark exchanges may receive contributions 
from new particles. As in the Standard Model E(m t ) « 0.3, the contribu- 
tion of Qq to e'/e is dominated by "11" which results from the operator 
mixing between Qq and other operators, in particular the current-current 
operator Q2. The latter mixing taking place at scales below Mw is unaf- 
fected by new physics contributions which can enter ( 4.19Q only through 



the function E. One would need an order of magnitude enhancement of E 
through new physics contributions in order to see a 40% effect in the con- 
tribution of Qq to e'/e and unless the sign of E is reversed one would find 
rather a suppression of e'/e with respect to the SM expectations than a re- 
quired enhancement. This discussion shows that a substantial enhancement 
of QCD-penguins through new physics contribution as suggested recently 
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in ||73[ is rather implausible. On the other hand in the case of the operator 



Qs one has instead of ([4.19 



(eye)- « [-10] Z(m u ..) (4.21) 



where _R 8 is given by (|4.20| ) with B^ 2 ^ replaced by B^ 2 ^ and Z(m t , ...) re- 



sults from electroweak penguin diagrams, in particular Z-penguin diagrams. 
Again "—10" comes from operator mixing under renormalization and new 
physics can enter Q4.21| ) only through the function Z. But this time any 
sizable impact of new physics on the function Z translates directly into a 
sizable impact on e'/e. 



5 Summary 

As we have seen, the estimates of e' / e in the Standard Model are typically 
below the experimental data. However, as the scanning analyses show, 
for suitably chosen parameters, e'/e in the Standard Model can be made 
consistent with data. Yet, this happens only if all relevant parameters 
are simultaneously close to their extreme values. This is clearly seen in 
table 0. Moreover, the probability density distributions for e'/e in fig. [1] 
indicates that values of e'/e in the ball park of the experimental grand 
average 21.2 ■ 10~ 4 are rather improbable. 

In spite of a possible "disagreement" of the Standard Model with the 
data, one should realize that certain features present in the Standard Model 
are confirmed by the experimental results. Indeed, the sign and the order 
of magnitude of e'/ e predicted by the SM turn out to agree with the data. 
In obtaining these results renormalization group evolution between scales 
C(M W ) and (9(1 GeV), an important ingredient in the evaluation of the 
Wilson coefficients Vi{n), plays a crucial role. As analyzed in [OJ without 
these renormalization group effects ye and y$ would be tiny and e'/e at 
most O(10~ 5 ) in vast disagreement with the data. Since these effects are 
present in all extensions of SM, we conclude that we are probably on the 
right track. 

Unfortunately, in view of very large hadronic and substantial parametric 
uncertainties, it is impossible to conclude at present whether new physics 
contributions are indeed required to fit the data. Yet as we stressed above, 
there is still a lot of room for non-standard contributions to e'/ e. The most 
plausible sizable new contributions could come from chromomagnetic pen- 
guins in general supersymmetric models and modified Z-penguins. On the 
other hand substantial modification of QCD-penguins through new physics 
are rather implausible. 

In view of large hadronic uncertainties it is difficult to conclude what 
is precisely the impact of the e'/e-data on the CKM matrix. However, 
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as analyzed in |]T3[ there are indications that the lower limit on lm\ t is 
improved. The same applies to the lower limits for the branching ratios for 
Kl — > ifivv and Ki — > 7r°e + e~ decays. 

The future of e'/ e in the Standard Model and in its extensions depends 
on the progress in the reduction of parametric and in particular hadronic 
uncertainties. In this context it is essential to get full control over renormal- 
ization scheme and renormalization scale dependences of hadronic matrix 
elements. Personally I believe that the only hope for making Standard 
Model naturally consistent with the data, without streching all parame- 
ters to their extreme values, is the value of m s (2 GeV) in the ball park of 
90 - 100 MeV, b£ /2) in the ball park of 1.5 - 2.0 and Bf /2) < 1.0 (both in 
the NDR scheme). The required values for m s seem to be found in lattice 
calculations but the story is not finished yet. Values of B^ 2 ^ in the ball 



park of 1.5 — 2.0 are suggested by the chiral quark model calculations [M 
and a higher order term 0(p 2 /N) in the large-N approach f54 |. Yet, from 



my point of view, it is not clear how well the chiral quark model approxi- 
mates QCD and whether other higher order terms in the large-N approach 
will weaken the indicated enhancement of B$ 2 \ We should hope that the 
new efforts by the lattice community will help in clarifying the situation. 
An interesting work in this direction in the framework of the large-N ap- 
proach can also be found in pEBfl . In any case e'/e already played a decisive 
role in establishing direct CP violation in nature and its rather large value 
gives additional strong motivation for searching for this phenomenon in 
other decays. 
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